Equilibrium Distributions in Open and Closed Statistical Systems 
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Abstract: In this communication, the derivation of the Boltzmann-Gibbs and the Maxwellian distributions is presented from 
a geometrical point of view under the hypothesis of equiprobability. It is shown that both distributions can be obtained by 
working out the properties of the volume or the surface of the respective geometries delimited in phase space by an additive 
constraint. That is, the asymptotic equilibrium distributions in the thermodynamic limit are independent of considering open 
or closed homogeneous statistical systems. 
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1 Introduction 

In this paper, different classical results 0] |2l are ob- 
tained from a geometrical interpretation of different 
multi-agent systems evolving in phase space under the 
hypothesis of equiprobability. 

We start by deriving in section [2] the Boltzmann- 
Gibbs (exponential) disuibution by means of the 
geometrical properties of the volume of an N- 
dimensional pyramid. The same result is obtained 
when the calculation is performed over the surface of 
a such ./V-dimensional body. In both cases, the moti- 
vation is a multi-agent economic system with an open 
or closed economy, respectively. 

Also, the Maxwellian (Gaussian) distribution is 
derived in section [3] from geometrical arguments over 
the volume or the surface of an iV-sphere. Here, the 
motivation is a multi-particle gas system in contact 
with a heat reservoir (non-isolated or open system) or 
with a fixed energy (isolated or closed system), re- 
spectively. 

Last section contains our conclusions. 



Derivation of the 
Gibbs Distribution 



Boltzmann- 



2.1 Multi-agent economic open systems 

Here we assume N agents, each one with coordinate 
Xi, i = 1, . . . , N, with Xi > representing the wealth 
or money of the agent i, and a total available amount 
of money E: 



Under random or deterministic evolution rules for the 
exchanging of money among agents, let us suppose 
that this system evolves in the interior of the N- 
dimensional pyramid given by Eq. £[]). The role of 
a heat reservoir, that in this model supplies money in- 
stead of energy, could be played by the state or by 
the bank system in western societies. The formula for 
the volume Vn(E) of an equilateral iV-dimensional 
pyramid formed by N + 1 vertices linked by N per- 
pendicular sides of length E is 



V N {E) 



E^_ 



(2) 



We suppose that each point on the iV-dimensional 
pyramid is equiprobable, then the probability 
f(xi)dxi of finding the agent i with money Xj is 
proportional to the volume formed by all the points 
into the (N — 1) -dimensional pyramid having the ith- 
coordinate equal to Xj. We show now that /(xj) is the 
Boltzmann factor (or the Maxwell-Boltzmann distri- 
bution), with the normalization condition 

f(xi)dxi = 1. (3) 

If the zth agent has coordinate Xj, the N — 1 re- 
maining agents share, at most, the money E — Xj on 
the (N — 1) -dimensional pyramid 

xi + x 2 h Xi-i + Xj + i h xn < E - Xi, (4) 

Cj). It can be easily 



whose volume is Vjy-i(E 
proved that 



X\ + X2~\ h Xat_i + Xn < E. 



(1) 



V N {E) = / V N -i(E - Xi )dxi 
Jo 



(5) 



Hence, the volume of the iV-dimensional pyramid 
for which the ith coordinate is between xi and Xj + 
dxi is Vn-\(E — xi)dxi. We normalize it to satisfy 
Eq. ©, and obtain 



f(xi 



V N -i(E-Xi) 



V N (E) ' 
whose final form, after some calculation is 

f{ Xl ) = NE-\l- X -±) N -\ 



(6) 



(7) 



If we call e the mean wealth per agent, E = iVe, then 
in the limit of large N we have 



lim ( 1 

iv>i 



X i 
E 



N-l 



-Xi/e 



(8) 



The Boltzmann factor e _Xi//e is found when N 3> 1 
but, even for small N, it can be a good approximation 
for agents with low wealth. After substituting Eq. ([8]) 
into Eq. ©, we obtain the Maxwell-Boltzmann distri- 
bution in the asymptotic regime N — > oo (which also 
implies E — > oo): 



f(x)dx = -e x ^dx, 



(9) 



where the index i has been removed because the distri- 
bution is the same for each agent, and thus the wealth 
distribution can be obtained by averaging over all the 
agents. This distribution has been found to fit the real 
distribution of incomes in western societies O. 

This means that the geometrical image of the 
volume-based statistical ensemble allows us to re- 
cover the same result than that obtained from the mi- 
crocanonical ensemble [4] that we show in the next 
section. 

2.2 Multi-agent economic closed systems 

Here, we derive the Boltzmann-Gibbs distribution by 
considering the system in isolation, that is, a closed 
economy. Without loss of generality, let us assume N 
interacting economic agents, each one with coordinate 
Xi, i = 1, . . . , N, with X{ > 0, and where Xj repre- 
sents an amount of money. If we suppose now that the 
total amount of money E is conserved, 



X\+ X 2 ~\ h XN-l + XN = E, 



(10) 



then this isolated system evolves on the positive part 
of an equilateral iV-hyperplane. The surface area 
Sn{E) of an equilateral iV-hyperplane of side E is 
given by 



S N (E) 



N 



(N-iy. 



E 



N-l 



(11) 



Different rules, deterministic or random, for the ex- 
change of money between agents can be given 
Depending on these rules, the system can visit the 
iV-hyperplane in an equiprobable manner or not. If 
the ergodic hypothesis is assumed, each point on the 
iV-hyperplane is equiprobable. Then the probability 
f(xi)dxi of finding agent i with money Xi is pro- 
portional to the surface area formed by all the points 
on the iV-hyperplane having the ith-coordinate equal 
to X{. We show that /(xj) is the Boltzmann factor 
(Boltzmann-Gibbs distribution), with the normaliza- 
tion condition ©. 

If the ith agent has coordinate Xj, the N — 1 re- 
maining agents share the money E—xi on the {N— 1)- 
hyperplane 

x\ + X2 h a?i-i + x i+ i h xn = E — Xi, (12) 

whose surface area is Sn-i(E — Xj). If we define the 
coordinate 9n as satisfying 



'JV-1 



sin on 



it can be easily shown that 



N 



S N (E)= f E S N ^(E- Xi )-^-, 
Jo sin 6 N 



(13) 



(14) 



Hence, the surface area of the iV-hyperplane for 
which the ith coordinate is between x,; and Xj + dxi 
is proportional to Sn-i(E — Xj)dxij sin 6 n- If we 
take into account the normalization condition Q, we 
obtain 



f(Xi 



1 S N -i(E-Xj) 
Sn(E) sin#/v 

whose form after some calculation is 

f{ Xl ) = {N-l)E-\l- X ^) N '\ 



(15) 



(16) 



If we call e the mean wealth per agent, E = Ne, then 
in the limit of large N we have 



lim 

AT>1 



1 



Xi 

E 



N-2 



-Xi/e 



(17) 



As in the former section, the Boltzmann factor e~ Xi ^ e 
is found when N 3> 1 but, even for small N, it can 
be a good approximation for agents with low wealth. 
After substituting Eq. (fTTT ) into Eq. ([TBI ), we obtain 
the Boltzmann distribution (O in the limit N — > 00 
(which also implies E — > 00). This asymptotic re- 
sult reproduces the distribution of real economic data 
||3l and also the results obtained in several models of 



economic agents with random exchange interactions 
0. 

Depending on the physical situation, the mean 
wealth per agent e takes different expressions and in- 
terpretations. For instance, we can calculate the de- 
pendence of e on the temperature, which in the micro- 
canonical ensemble is defined by the derivative of the 
entropy with respect to the energy. The entropy can be 
written as S = —kN J °° f(x) In f(x) dx, where f{x) 
is given by Eq. Q and k is Boltzmann's constant. If 
we recall that e = E/N, we obtain 



E 

S{E) = kNln— + kN. 
The calculation of the temperature T gives 



dS_ 
~dE 



N 



kN _k 
~E~~~e' 



(18) 



(19) 



Thus e = kT, and the Boltzmann distribution is ob- 
tained in its usual form: 



f(x)dx 



1 

kf 



c ' kT dx. 



(20) 



3 Derivation of the Maxwellian Dis- 
tribution 

3.1 Multi-particle open systems 

Let us suppose a one-dimensional ideal gas of N non- 
identical classical particles with masses wij, with i = 
1, . . . , N, and total maximum energy E. If particle i 
has a momentum rriiVi, we define a kinetic energy: 



Pi = 2 m i V i> 



(21) 



where pi is the square root of the kinetic energy Ki. If 
the total maximum energy is defined as E = R 2 , we 
have 



Pi 



+ P 2 2 + ---+P 2 n~i+P 2 n<R 2 



(22) 



We see that the system has accessible states with dif- 
ferent energy, which can be supplied by a heat reser- 
voir. These states are all those enclosed into the vol- 
ume of the iV-sphere given by Eq. (1221 . The formula 
for the volume Vn (R) of an iV-sphere of radius R is 



V N (R) 



7T 2 



r(f + i) 



-R 



N 



(23) 



where T(-) is the gamma function. If we suppose that 
each point into the iV-sphere is equiprobable, then the 



probability fipijdpi of finding the particle i with co- 
ordinate pi (energy pf) is proportional to the volume 
formed by all the points on the iV-sphere having the 
fth-coordinate equal to p^. We proceed to show that 
f(pi) is the Maxwellian distribution, with the normal- 
ization condition 



f R 

/ f{Pi)dpi = 1. 
J-R 



(24) 



If the ith particle has coordinate p it the (N — 1) 
remaining particles share an energy less than the max- 
imum energy R 2 — p\ on the (N — l)-sphere 

Pi +Pl ■ ■ ■ +Pi-l +Pi+i ■ ■ ■ +Pn < R 2 -Pi (25) 



whose volume is Vjv-i(y R 2 ~ Pi)- ^ can be easily 
proved that 

V N (R) = [ R V N - l {^R 2 -p 2 )dp l . (26) 

Hence, the volume of the iV-sphere for which 
the ith coordinate is between pi and pi + dpi is 

Vn-i{\] R 2 — Pi)dpi- We normalize it to satisfy 
Eq. (|24l . and obtain 



V N . l UR 2 -p 2 ) 



V N (R) 

whose final form, after some calculation is 

f(p t ) = C N R- 1 (l-^) 2 , 



with 



Cn 



i r(gg) 



(27) 



(28) 



(29) 



For N S> 1, Stirling's approximation can be applied 
to Eq. d29l . leading to 



lim Cn ~ -h=\\ —■ 

tv»i V^V 2 



(30) 



If we call e the mean energy per particle, E = R 2 
Ne, then in the limit of large N we have 



JV-l 

,2 \ 



Af>l \ R 2 



(31) 



The factor e~ p ^ 2e is found when N S> 1 but, even for 
small N, it can be a good approximation for particles 
with low energies. After substituting Eqs. (l30l-(f3TT) 
into Eq. (|28T ). we obtain the Maxwellian distribution 



in the asymptotic regime N 

E -> oo): 



oo (which also implies 



then 



1 



f(p)dp=\l^- e e- p /2e dp, 



(32) 



where the index i has been removed because the dis- 
tribution is the same for each particle, and thus the ve- 
locity distribution can be obtained by averaging over 
all the particles. 

This newly shows that the geometrical image of 
the volume-based statistical ensemble allows us to re- 
cover the same result than that obtained from the mi- 
crocanonical ensemble [5] that it is presented in the 
next section. 

3.2 Multi-particle closed systems 

We start by assuming a one-dimensional ideal gas of 
iV non-identical classical particles with masses rrii, 
with i = 1, . . . , N, and total energy E. If particle i 
has a momentum m^i, newly we define a kinetic en- 
ergy Ki given by Eq. ([2TT> . where pi is the square root 
of Ki. If the total energy is defined as E = R 2 , we 
have 



2 i 2 i 
Pl+P2 + 



+ P%-1 +P% 



R 2 



(33) 



Rd9 



dpi 



(37) 



It can be easily proved that 

S N (R)= S N -i(R cos 6) Rd6. (38) 

J-n/2 

Hence, the surface area of the iV-sphere for which 
the ith coordinate is between and pi + dpi is 
SjV-iX-Rcos 9)Rd9. We rewrite the surface area as 
a function of pi, normalize it to satisfy Eq. (f24l> . and 
obtain 



f{Pi) 



1 S N ^R 2 -p 2 
S N {R) (1 _ 4)V2 

^ Ff. ' 



(39) 



whose final form, after some calculation is 

f(p l ) = C N R- 1 (l-^) 2 , 



with 



Cn 



1 r(f) 
^r(^i)- 



(40) 



(41) 



For N 3> 1, Stirling's approximation can be applied 
to Eq. (HTT ). leading to 



We see that the isolated system evolves on the sur- 
face of an iV-sphere. The formula for the surface area 
Sn{R) of an iV-sphere of radius R is 



Sn(R) 



N_ 

2vr 2 rN 



r(f) 



(34) 



where T(-) is the gamma function. If the ergodic 
hypothesis is assumed, that is, each point on the iV- 
sphere is equiprobable, then the probability f(pi)dpi 
of finding the particle % with coordinate pi (energy 
p 2 ) is proportional to the surface area formed by all 
the points on the iV-sphere having the ith-coordinate 
equal to p^. Our objective is to show that f(pi) is the 
Maxwellian distribution, with the normalization con- 
dition (1241 . 

If the ith particle has coordinate pi, the (N — 1) 
remaining particles share the energy R 2 — p 2 on the 
(N — l)-sphere 



Pi +Pl ■ ■ ■ +P 2 i-i +P 2 i+i ■ ■ ■ +Pn 



R 2 



■Pi, 



(35) 



whose surface area is Sjv-i(y R 2 — P 2 )- If we define 
the coordinate 6 as satisfying 



R 2 cos 2 9 = R 2 



ph 



(36) 



lim Cn ~ -^=\ —. 

iv»i a/ttV 2 



(42) 



If we call e the mean energy per particle, E = R 2 
Ne, then in the limit of large N we have 



N-3 
J2 \ ~^2~ 



lim 1 - ^ 

w>i \ R 2 



(43) 



As in the former section, the Boltzmann factor e _p * / 2e 
is found when N 3> 1 but, even for small N, it can be 
a good approximation for particles with low energies. 
After substituting Eqs. ([42]>-([43]> into Eq. (00]), we ob- 
tain the Maxwellian distribution (1321 in the asymp- 
totic regime N — > 00 (which also implies E — > 00). 

Depending on the physical situation the mean en- 
ergy per particle e takes different expressions. For an 
isolated one-dimensional gas we can calculate the de- 
pendence of e on the temperature, which in the micro- 
canonical ensemble is defined by differentiating the 
entropy with respect to the energy. The entropy can 
be written as S = —kN f(p) In f(p) dp, where 
f(p) is given by Eq. (1321 and k is the Boltzmann con- 
stant. If we recall that e = E/N, we obtain 

S (E) = hN In (^\+hN(ln(2ir)-l). (44) 



The calculation of the temperature T gives 



T 



-i 



dS\ kN _ k 
dEJ N ~ 1 2E ~2~e' 



(45) 



Thus e = kT/2, consistent with the equipartition the- 
orem. If p 2 is replaced by \mv 2 , the Maxwellian dis- 
tribution is a function of particle velocity, as it is usu- 
ally given in the literature: 



9(v)dv = J-^-e-™ 2 / 2kT dv. (46) 



2nkT 



4 Conclusion 

We have shown that the Boltzmann factor describes 
the general statistical behavior of each small part of a 
multi-component system whose components or parts 
are given by a set of random variables that satisfy an 
additive constraint, in the form of a conservation law 
(closed systems) or in the form of an upper limit (open 
systems). The derivation of this factor for open sys- 
tems in a general context has been presented in |6). 

Let us remark that these calculations do not need 
the knowledge of the exact or microscopic random- 
ization mechanisms of the multi-agent system in or- 
der to reach the equiprobability. In some cases, it can 
be reached by random forces O, in other cases by 
chaotic Q or deterministic [8] causes. Evidently, the 
proof that these mechanisms generate equiprobability 
is not a trivial task and it remains as a typical chal- 
lenge in this kind of problems. 

In summary, this work has presented a straightfor- 
ward geometrical argument that recalls us the equiv- 
alence between canonical and microcanonical ensem- 
bles in the thermodynamic limit in the particular con- 
text of physical sciences. For the general context of 
homogeneous multi-agent systems, we conclude by 
highlighting the statistical equivalence of the volume- 
based and surface-based calculations in this type of 
systems. 
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